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Abstract 

We present a study on physical properties for a large distant red galaxy 
(DRG) sample, using the K-selected multi-band photometry catalog of the 
GOSMOS/UltraVISTA held and the GANDELS NIR data. Our sample includes 4485 
DRGs with (J — iF)AB > 1-16 and Aab <23.4 mag, and 132 DRGs have HST/WFG3 
morphological measurements. The results of nonparametric measurements of DRG 
morphology are consistent with our rest-frame UVJ color classihcation: quiescent 
DRGs are generally compact while star-forming DRGs tend to have extended struc¬ 
tures. We hnd the star formation rate (SFR) and the stellar mass of star-forming 
DRGs present tight “main sequence” relations in all redshift bins. Moreover, the spe- 
cihc SFR (sSFR) of DRGs increase with redshift in all stellar mass bins and DRGs 
with higher stellar masses generally have lower sSFRs, which indicates that galaxies 
were much more active on average in the past, and star formation contributes more 
to the mass growth of low-mass galaxies than to high-mass galaxies. The infrared 
(IR) derived SFR dominate the total SFR of DRGs which occupy the high-mass 
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range, implying that the J — K color criterion effectively selects massive and dusty 
galaxies. DRGs with higher M* generally have redder (G —l/)rest colors, and the 
{U — ld)rest colors of DRGs become bluer at higher redshifts, suggesting high-mass 
galaxies have higher internal dust extinctions or older stellar ages and they evolve 
with time. Finally, we find that DRGs have different overlaps with EROs, BzKs, 
lEROs and high-z ULIRGs indicating DRGs is not a special population and they can 
also be selected by other color criteria. 

Key words: galaxies: evolution - galaxies: formation - galaxies: high-redshift - 
galaxies: structure 


1. Introduction 

The universe at high redshift looks very different from that at low redshift and local. At 
the epoch of 2 ; ~ 2, a large fraction of massive galaxies are still active, and their morphology are 
quite different from that of nearby galaxies (McGracken et ah 2010; Fang et al. 2012; Wang et ah 
2012). In order to examine current theoretical models of galaxy formation and evolution, studies 
of galaxies over a wide range of lookback time are needed (Shapley 2011). Within the past 
decade, many novel techniques have been applied to assemble multi-wavelength observations of 
galaxies at this important epoch and significant progress has been made in our understanding 
of high-redshift galaxies (Steidel et al. 1996; Small et al. 1997; Franx et al. 2003; Daddi et al. 
2004; Gronwall et al. 2007; Yan et al. 2007; Dey et al. 2008; Dunne et al. 2009; Gaputi et al. 
2012; Arcila-Osejo & Sawicki 2013; Ilbert et al. 2013). 

At early time, one method for selecting high-redshift galaxies is the Lyman break tech¬ 
nique, which uses the significant break at 912 A in the galaxy rest-frame UV spectrum to 
select star-forming galaxies (so-called Lyman Break Galaxies, or LBGs), and was first applied 
at 2 ; ~ 3 (Steidel et al. 1996). The LBGs are UV-bright star-forming galaxies and dominated by 
relatively low-mass systems (~ 10^*^ Mq). Gonsider that most of local elliptical and SO galax¬ 
ies commonly have stellar masses great than IO^^Mq, the LBGs may not be the progenitors 
of these massive galaxies today (Gonselice et al. 2007), and the Lyman break technique have 
missed a significant population of red massive galaxies at high redshift, such as distant red 
galaxies (DRGs) (Franx et al. 2003), BzKs (Daddi et al. 2004), ultraluminous infrared galaxies 
(ULIRGs) (Ghapman et al. 2003), etc. As a supplement to these UV-bright LBGs, Franx et 
al. 2003 introduced a new photometric technique for selecting z>2 massive red galaxies. The 
relevant spectral break is the Balmer break at 3650 A which attributes to the spectral shape of 
young stars and the 4000 A break which is caused by the absorption of metals in old stars. Since 
these breaks shift into the J band at 2 ; ~ 2, a red J — K color is a simple and effective criterion 
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(Franx et al. 2003). They simulated the tracks in the J — K color versus redshift diagram with 
model spectra taken from Bruzual & Chariot 1993, and imposed a criterion (J —iF)vega > 2.3 
for selecting galaxies at z > 2 (so-called Distant Red Galaxies, or DRGs). 

DRGs commonly have higher star formation rates (SFRs) (~ lOOMQ/yr) and more dust 
extinctions {Ay = 1 — 3 mag) (Franx et al. 2003; Forster Schreiber et al. 2004; Papovich et 
al. 2005; Webb et al. 2006). The redshift distribution of DRGs is quite broad, and the (J — 
-^)vega > 2.3 selection criterion does not uniquely sample high-mass galaxies (Papovich et al. 
2005; Grazian et al. 2006; Conselice et al. 2007). However, DRG selection technique does locate 
very massive galaxies ( 10^^ at z > 2 (van Dokkum et al. 2006; Conselice et al. 2007). 

In the present-day Universe, most of massive (M* > IO^^Mq) galaxies are quiescent systems. 
On the contrary, at 2 ; ~ 2 a large proportion of massive galaxies are still active, while others are 
“red and dead” with little or no ongoing star-forming activities (Labbe et al. 2005; Papovich et 
al. 2006; Daddi et al. 2007; Wang et al. 2012). Transformation from the active systems into the 
evolved systems is a reasonable explanation to this dramatic evolution. Some studies suggested 
that DRGs show properties of both passive systems and dust obscured active systems (Forster 
Schreiber et al. 2004; Labbe et al. 2005; Papovich et al. 2005; Webb et al. 2006). Labbe et 
al. 2005 showed color tracks for stellar population models (Bruzual & Chariot 2003) in the 
{I — K) vs. {K — [4.5]) diagram, and found that 70% of DRGs are well £t with dusty star¬ 
forming models and 30% are best described by old stellar population models. Webb et al. 
2006 combined Spitzer 24/im flux, and found a similar percentage (~65%) of DRGs are dusty 
active galaxies. By studying the spectra, Kriek et al. 2006a; Kriek et al. 2006b found that 
the stellar populations differ greatly among DRGs, from dusty starburst with small break to 
evolved galaxy with strong break and no line emission, and both the Ha measurements and 
stellar continuum modeling imply 45% of iC-selected massive galaxies at 2 ~ 2 are not forming 
stars intensely. 

The morphologies of DRGs are diverse. Toft et al. 2005 analyzed the morphologies of 
DRGs in the Hubble Ultra Deep Field based on HST Advanced Camera for Surveys (ACS) 
images and Near Infrared Camera and Multi-Object Spectrometer (NICMOS) images, and 
found that the rest-frame optical morphology of DRG is quite different from the rest-frame UV 
morphology. In the near-infrared images, some DRGs are extended and others are compact. 
Conselice et al. 2007 examined the morphologies of DRGs in the Extended Groth Strip (EGS) 
held based on ACS F814W images with eye-ball and nonparametric classihcation method, and 
found that the morphologies of DRGs show a diversity of types at all redshifts, from elliptical, 
peculiar to disk. 

The (J — A")vega > 2.3 criteria is a simple but outstanding selection method utilising 
the 3650 and 4000 Abreaks for high- 2 ; massive red galaxies, which are supplements to those 
UV-bright low-mass LBGs. The DRG sample provides us with the chance to investigate the 
nature of the population of the massive red objects at high redshifts, which are supposed to be 
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the progenitors of local massive populations, thus it makes us better understand the evolution 
of massive galaxies. Since the DRG samples of previous research programs are commonly 
small (< 1000), it still remains many uncertainties on the physical properties of distant red 
galaxies. In order to have a thorough understanding of this kind of galaxies, in this paper we 
present the largest sample of DRGs so far, which is selected from the COSMOS/UltraVISTA 
held, based on the high quality NIR data presented by McCracken et ah 2012 and Muzzin 
et ah 2013a. To clearly describe the details of the characteristics of DRGs, we investigate 
the correlations between their SFRs, redshifts, masses and morphologies. The rest-frame UV 
emission of galaxies is mainly contributed by the hottest stars and can be severely affected 
by dust extinction, therefore it is essential to study z ^ 2 galaxies from observed NIR bands, 
which probe the rest-frame optical morphologies. Our morphology analysis was enabled by 
the HST Wide Field Camera 3 (WFC3) NIR imaging of the Cosmic Assembly Near-infrared 
Deep Extragalactic Legacy Survey (CANDELS)(Grogin et al. 2011; Koekemoer et ah 2011), 
and we compare the nonparametric measure of galaxy morphology with visual inspection and 
UVJ color-color diagram. In the meanwhile, we also compare DRG with other color selection 
criteria. 

The paper is organized as follows. We introduce the iF-selected catalog of the 
COSMOS/UltraVISTA held in Section 2, and the DRG sample selection method and clas- 
sihcation in Section 3. We show our results on physical properties of DRGs in Section 4. We 
compare the DRG method with other color selection methods in Section 5, and conclude our 
results in Section 6. Throughout this paper, we assume an Dm = 0.3, Da = 0.7 and Hq = 70 
km s“^ Mpc“^ cosmology, and all magnitudes and colors are given in AB system unless stated 
otherwise. Symbol such as “[3.6]” means the AB magnitude at the wavelength 3.6 /im. 

2. Observations and data 

The COSMOS held is centered at a = 10'*00”^28.6", S = +02'^12™21* (J2000.0), and 
covers an area of ~ 2 square degree (Scoville et al. 2007). The COSMOS project gathers a 
variety of ground- and space-based observations, covering nearly the whole spectrum range 
(Martin et al. 2005; Capak et al. 2007; Basinger et al. 2007; Sanders et al. 2007; Schinnerer et 
al. 2007; Taniguchi et al. 2007; Lilly et al. 2009; Ilbert et al. 2010; McCracken et al. 2010). 

The multi-band photometry data we used in this paper come from the iF-selected catalog 
of the COSMOS/UltraVISTA held provided by Muzzin et al. 2013a, which is produced based 
on the near infrared (NIR) data from the UltraVISTA data release 1 (McCracken et al. 2012). 
The high-quality and deep-enough NIR data make it possible for us to select a large sample 
of DRGs. The catalog contains PSF-matched photometry in 30 photometric bands covering 
the wavelength from 0.15/im to 24/im, and includes the available GALEX (Martin et al. 2005), 
CFHT/Subaru (Capak et al. 2007), UltraVISTA (McCracken et ah 2012), S-COSMOS (Sanders 
et al. 2007), and zCOSMOS (Lilly et al. 2009) datasets. After removing regions contaminated 
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by bright stars, the effective area of overlap between optical and NIR data is 1.62 deg^. The 
photometry in all bands is corrected for Galactic dust attenuation, and stars are separated 
from galaxies using the (J — K) vs. {u — J) color space, and it agrees well with SExtractor’s 
class_star parameters. The catalog is presented as a set of fluxes in the 2"1 color aperture 
with an AB zeropoint of 25.0, and the 90% completeness limit of the catalog is A"s,tot = 23.4. 

The derived physical parameters we used in this paper also come from the catalog 
provided by Muzzin et al. 2013a. Photometric redshifts are derived for all galaxies using the 
EAZY code (Brammer et al. 2008) with 30 bands, and have been calibrated by the highest-quality 
spectroscopic redshifts from zCOSMOS. The stellar population parameters such as stellar mass 
(M*), age, Ay, SFRsed are derived by SED htting method in FAST code (Kriek et al. 2009a). 
The rest-frame UV luminosity is derived by EAZY. The total IR luminosity is derived by 
fitting the 24/im flux using the log-average of the Dale & Helou 2002 templates, 
but note that using the log-average of the Chary & Elbaz 2001 templates would 
provide very similar results. Then the rest-frame UV and the total IR luminosities are 
converted to SFRuv,uncorr and SFRir according to Kennicutt 1998. The total star formation 
rate of the galaxy can then be determined via SFRtot = SFRuv,uncorr+ SFRir. The rest-frame 
U — V and V — J colors are calculated for all galaxies by EAZY, which determines the colors 
by integrating the best-ht SED through the redshifted Alter curves over the appropriate wave 
length range. 

3. Selection and classification of DRGs 
3.1. DRG Sample Selection 

To select a sample of DRGs from Ultra VISTA catalog, flrstly we remove all the objects 
classifled as stars and contaminated galaxies deflned by Muzzin et al. 2013a. In order to make 
our photometric selecting criterion close to the original DRG selecting criterion introduced by 
Franx et al. 2003, we make a simulation in BG03 code (Bruzual & Gharlot 2003) to determine 
the corrections for J and K bands. Finally, we set a small correction for J — K color: (J — 
A')uitraVISTA + 0.14 = (J — IF)ISAAC (AB), SO We set (J — iF)uitraVISTA > 1-16 (AB), which is 
equivalent to (T —iF)isAAc > 1-3 (AB), as the criterion for selecting DRGs. On the other hand, 
we carefully compared the location of the stellar sequence with previous works, and And the 
distance on J — iF color between stellar sequence and DRG selecting criterion is 1.36, which is 
consistent with previous works (Foucaud et al. 2007; Kajisawa et al. 2008). 

Figure 1 shows the JF-selected galaxies (objects in masked area are removed) in the 
[J — K) vs. K diagram, the dashed line represents J — K = 1.16, which is our DRG selecting 
criterion, objects above which are all DRGs and are plotted in red color. The dot-dashed line 
represents the stellar sequence, and all stars are plotted in green color. The total number of 
DRGs we select from the GOSMOS/UltraVISTA held is 4485 with K < 23.4 and J < 24.2 mag 
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Fig. 1. [J — K) color vs. K magnitude for galaxies in the COSMOS/UltraVISTA field. The dashed line 
corresponds to J —A" = 1.16 which is our DRG selecting criterion, galaxies above which are all DRGs 
and are plotted in red color. The dot-dashed line represents the stellar sequence, and stars are plotted in 
green. The dot-dot-dashed line indicates the depth of the J band. 

(Therefore, all DRGs in our sample actually have K < 24.2 — 1.16 = 23.04 mag). 

Figure 2 shows the iF-band differential number counts of DRGs and all iF-selected galax¬ 
ies. For comparison, we have also plotted the number count results from other literature. In 
this figure, the red squares and the black triangles represent the number counts of DRGs and all 
R-selected galaxies in the GOSMOS/UltraVISTA field, respectively. The blue stars represent 
the number counts of DRGs whose J-band magnitude limit is not considered. Results of DRGs 
from GOODS-S (Grazian et ah 2006), HDF-S (Labbe et al. 2003), GOODS-N (Kajisawa et al. 
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Fig. 2. K-h&nd differential number counts for DRGs (red squares) and all iiT-selected galaxies (black 
triangles) in the COSMOS/UltraVISTA field. The blue stars represent the number counts of DRGs 
regardless of their bu detection in J-band or lack there of. For comparison, we also overplot the number 
counts of DRGs and all isT-selected galaxies from previous works. 
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2006), UDS-EDR & DXS-EDR (Foucaud et al. 2007), GOODS-N (Reddy et al. 2005), UKIDSS 
(Kim et al. 2011) and EGS (Barro et al. 2011) are also plotted. And results of all JGselected 
galaxies from Deep3a-F & Daddi-F (Kong et al. 2006), GOSMOS (McGracken et al. 2010) and 
UKIDSS (Kim et al. 2011) are plotted in different symbols as well. As shown by this hgure, 
our number counts of DRGs and all Tr-selected galaxies in the GOSMOS/UltraVISTA held are 
in good agreement with those in previous works in other helds. 

3.2. Classification based on rest-frame UVJ colors 

As shown by previous literatures (Labbe et al. 2005; Webb et al. 2006; Kriek et al. 
2006a; Kriek et al. 2006b), DRG is a diverse population, most of the DRGs are dusty star¬ 
forming galaxies while others are old and red with little on going star formation activities. In 
order to specify the details of the physical properties and the evolution of them, separating 
DRGs into star-forming and quiescent populations is necessary. 

We classify our sample as star-forming DRGs (sDRGs) and quiescent DRGs (qDRGs) 
based on the rest-frame UVJ color criteria dehned by Muzzin et al. 2013b. According to their 
dehnition, qDRGs are classihed by: 

U — V>1.3, U — J< 1.5, [all redshifts], (1) 

U-V>{V-J)x 0.88 + 0.69, [0.0 < < 1.0], (2) 

U-V>{V-J)x 0.88 + 0.59, [1.0 < z < 4.0]. (3) 

Otherwise, they are star-forming DRGs. These classihcation criteria are plotted as dashed 
lines in Figure 3. Based on this classihcation method, we identify 3725 sDRGs and 760 qDRGs, 
the percentage of sDRGs is 83%, which is a little higher than previous results (Labbe et al. 
2005; Webb et al. 2006). 

3.3. Classification based on visual morphologies 

Morphologies of DRGs provide direct information on the formation and evolution history 
of these objects, and correlate a range of physical properties such as stellar mass, SFR and 
rest-frame color. It is difficult to study morphologies of high redshift dusty galaxies based 
on their observed optical images, because their observed optical light probes the rest-frame 
UV emission for objects at z ~ 2, and their apparent morphologies can be easily changed by 
patchy dust extinction. For example, Gonselice et al. 2007 showed that DRGs have peculiar, 
clumpy morphologies in the HST /AGS F814W band image. Based on the HST NIGMOS-Ji7 
(0."09/pixel) and AC^-BViz (0."03/pixel) images, however. Toft et al. 2005 found that the rest- 
frame optical morphology of DRGs is quite different from the rest-frame UV morphology. The 
rest-frame UV emission of galaxies mainly contributed by the hottest stars and can be severely 
affected by dust extinction, therefore it is essential to study z ^2 galaxies from observed NIR 
bands which probe the rest-frame optical morphologies. Our morphology analysis was enabled 
by the i7S'T/WFG3 NIR imaging with a resolution of 0''06/pixel from the Gosmic Assembly 
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Fig. 3. Classification for DRGs based on rest-frame {U — V) vs. {V — J) colors in six redshift bins. This 
color criteria (dashed lines) used to separate star forming DRGs (sDRGs) and quiescent DRGs (qDRGs) 
are dehned by Muzzin et al. 2013b. The number of DRGs, the fraction of sDRGs and qDRGs in each 
redshift bins are shown, respectively. 


Near-infrared Deep Extragalactic Legacy Survey (CANDELS) (Grogin et al. 2011; Koekemoer 
et al. 2011). 

Since the redshift distribution of DRGs is quite broad, we analyze their rest-frame optical 
morphologies on WFG3 F125W (J) or F160W {H) bands according to their redshifts. For 
DRGs with 1 < z < 1.6, we choose WFG3 F125W (J) bandpass for morphological analysis, it 
corresponds approximately to R-band in the rest-frame in this redshift range, but in the redshift 
range of 1.6 < ^ < 3, we analyze galaxy morphology in the rest-frame V from the F160W (II) 
image instead. Finally, 37 DRGs in our sample have J-band counterparts (1 < z < 1.6), and 
95 DRGs (1.6 < 2 ; < 3) are detected at Lf-band image. The sensitivity (Sa) of the GANDELS 
imaging data reach 27 mag in the F125W and the F160W, and all of DRGs in our sample are 
brighter than 24 mag in both J and H bands, thus the GANDELS data is deep enough to 
morphologically classify this DRG sample. 

We perform the visual inspection for this sample by three of us independently, we com¬ 
bined each classiher’s results and review the images together to resolve the disagreements. The 
disagreements cover about 10% of this sample, and mainly exist among galaxies which can’t 
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Table 1. Numbers of star-forming and quiescent DRGs, classified by rest-frame UVJ colors in three visual morphological 
types. 


DRGs (1 < 0 < 3) 

Spheroid 

Disk 

Irregular 

Total 

Star-forming DRGs (sDRGs) 

20 

33 

48 

101 

Quiescent DRGs (qDRGs) 

24 

5 

2 

31 


Spheroids Disks Irregulars 



Fig. 4. Examples of three general morphological types from visual inspection: spheroid on the 
left, disk in the middle and irregular on the right panel. The size of each image is 2'.'4 x 2f.'4. 

be clearly classified as disks or irregular galaxies. We find that the DRGs in our sample show 
very diversified morphologies, covering a wide range of types from extended disks, compact 
spheroids to clumpy and irregular morphologies. We define three general morphological types 
as follows: 

Spheroid : round and centrally concentrated source with no extended structure. 

Disk : single bright core with smooth extended disk-like structure. 

Irregular : clumpy structure with more than two bright cores, sometimes shows interaction 
features such as tidal arms. 

Examples of three morphological types are shown in Figure 4. The result of the visual 
inspection supports the finding that the significant evolution of the Hubble sequence has already 
existed at high redshift (Kajisawa & Yamada 2001). In Table 1 we list the numbers of star¬ 
forming and quiescent DRGs classified by rest-frame UVJ colors in three visual morphological 
types. Among star-forming DRGs, the spheroids occupy a proportion of 19.8%, disks occupy 
32.7% and irregulars make up 47.5%. And among quiescent DRGs, the spheroids make up 
77.4%, disks make up 16.1% and irregulars make up 6.5%. We find that most of quiescent 
DRGs are compact and star-forming DRGs tend to have extended structures. 
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Fig. 5. Photometric redshift distribution (left panel) and stellar mass distribution (right panel) of DRGs. 

All DRGs are plotted in black lines, while star-forming and quiescent DRGs classified by rest-frame UVJ 
color are plotted in blue and red, respectively. 

4. Results 

4 . 1 . Redshift and mass distributions 

The photometric redshifts from the Ultra VISTA catalog are in good agreement with the 
spectroscopic redshifts from zCOSMOS (Lilly et ah 2009), with an average dz/{l + z) = 0.013 
(Muzzin et ah 2013a). We show the redshift distributions of both star-forming and quiescent 
DRGs in the left panel of Figure 5. The star-forming DRGs (blue lines) distribute at 0 < z < 3.5 
with {z) = 1.74, and the quiescent DRGs (red lines) are in the range 2 ; ~ 1 — 3 with {z) = 1.98. 
For all DRGs, our sample has {z)=1.79. In total, there are 3% DRGs distribute at z <1, 95.5% 
at 1 < z < 3 and 1.5% at z > 3. The right panel of Figure 5 shows the stellar mass distributions 
of both star-forming and quiescent DRGs. We hnd the masses of them are relatively high, and 
mainly range from IO^^Mq to 10 ^^'^Mq. The mean stellar masses are 10 ^°''^Mq, 10^°'®Mq and 
10^° ®Mq for all, star-forming and quiescent DRGs in our sample, respectively. 

We compare the redshift distributions of DRGs with those of all of galaxies in four 
stellar mass bins in Figure 6. We hnd that in higher mass bins (e.g., IO^^'^Mq — IO^^Mq, 
IO^^Mq — IO^^Mq), DRGs represent nearly all of galaxies at z > 2. At M* > IO^^'^Mq and 
z > 2, the fraction of DRGs out of all galaxies is 70%, which is in good agreement with van 
Dokkum et ah 2006. From the analysis above, we demonstrate that although DRG selecting 
criterion also sample low-mass galaxies at low redshifts, DRGs occupy large proportion of 
galaxies at high-mass and high- 2 ; range, implying that DRG selecting criterion is sensitive to 
massive galaxies at high redshifts. Our results are also demonstrated in Gonselice et al. 2007. 
The spectroscopy study of Kriek et al. 2008 also supports our hnding that red galaxies dominate 


11 


































Num. 



c) 10.5<logio(M*/Mp,)<11.0 d) 11.0<logio(M*/MQ)<12.0 

I I All galaxies. 
\///\ non—DRGs 
^ DRGs 



^ph 


^ph 


Fig. 6. Redshift distributions of all galaxies (black lines), DRGs (red lines) and non-DRGs (green 
lines) in the COSMOS/UltraVISTA field in four stellar mass bins. 


12 


200 400 600 0 200 400 600 












































































' 1 ' r 

■ a) F125W 1.0<z<1.6 



1 1 1 1 1 

■ b) F160W 1.6<z<3.0 



sDRGs 


sDRGs 


A 

qDRGs 


qDRGs 

1 

* ‘ A 


1 

“A “ 


A A A 



* A 


★ 



A 

0 

* A* 


0 

A 

A A ^ A 

CM 

s 

\ ★ ★ 


CM 

s 

► 

► 

► 

► 

♦ 

♦ 

in 

1 

A * * 

A**‘ ** 

♦ 

IT) 

1 

A ♦ ♦♦♦ 

- * ♦ ♦ ♦ ♦ ■ 


A 




★ * 

A Irregular ★ 

♦ 


* *A *** * 

A Irregular a * ^ ♦ 

CM 

_ * Disk * 


CM 

★ Disk 


. .* Spheroid 

_ 1 _ 1 _ 1 _L 



. * Spheroid 

_ 1 _ 1 _ 1 _ 1 _ 1 _ 


0.4 0.6 0.4 0.6 

Gini Gini 


Fig. 7. M 20 vs. Gini coefficient for DRGs. The star-forming DRGs (sDRGs) are plotted in blue 
color, while quiescent DRGs are in red. Spheroid, disk and irregular DRGs which are classified by vi¬ 
sual inspection, are plotted in diamond, star and triangle symbols, respectively. DRGs on left panel 
with redshift 1 < z < 1.6, whose morphological parameters are measured from F125W image. DRGs on 
right panel distribute at 1.6 < z < 3, whose morphological parameters are measured from F160W image. 

the high-mass end of the galaxy population at z ~ 2 — 3. 

We have checked the mass completeness for our DRG sample, and hnd nearly all DRGs 
are above and the majority of them are above at 1 < z < 3. The distance 

between the majority of DRGs and the 90% completeness boundary at 1 < 2 ; < 3 are safe enough, 
which ensures our DRG is a complete sample. Only a very small proportion of DRGs are close 
to the completeness boundary, and most of them are located in low-mass and \ow-z range. 

4-2. Nonparametric morphology 

To further describe the morphologies of these sources, we have performed nonparametric 
measures of galaxy morphology, such as Gini coefficient (the relative distribution of the galaxy 
pixel flux values, or Gini) and M 20 (the second-order moment of the brightest 20% of the 
galaxy’s flux) (Lotz et ah 2006), using the Morpheus-software developed by Bob Abraham 
(Abraham et al. 2003). Figure 7 shows the distribution of DRGs on the Gini—M 20 plane in 
F125W and F160W bands, which corresponds to the redshift range 1 <z < 1.6 and 1.6 < z <3, 
respectively. Diamonds, stars and triangles represent spheroid, disk and irregular galaxies, 
classihed by visual inspection. Star-forming and quiescent DRGs classihed by the rest-frame 
UVJ colors are plotted in blue and red, respectively. We hnd that most of quiescent DRGs are 
visually compact, and have larger Gini and smaller M 20 values, while the star-forming DRGs 
contain different kinds of morphological types but most of them are visually extended and have 
smaller Gini and larger M 20 values. Similar results are also reported in Lee et al. 2013 and 
Bassett et al. 2013. 
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Kriek et al. 2009a; Kriek et al. 2009b support our morphological results that quiescent 
galaxies are commonly found to be ultra-compact, while galaxies with high SFRs have irregular 
and clumpy morphologies. Our morphological results also support the evolutionary tracks 
described in Barro et al. 2013, Bedregal et al. 2013 and Tadaki et al. 2014 that compact 
DRGs (with larger Gini and smaller M 20 ) contain compact quiescent galaxies (cQGs) and their 
progenitors compact star-forming galaxies (cSFGs). In our sample, there are 54.5% (24/44) 
spheroidal DRGs are quenched, namely they are cQGs, and the rest of them are cSFGs. The 
cSFGs are formed by gas-rich processes such as major mergers at high redshifts, and rapidly 
quenched (probably by AGN) into cQGs, most of which are classihed as spheroidal quiescent 
DRGs. These compact QGs at 2 ; ~ 2 may be the cores of today’s massive elliptical galaxies, as 
described by van de Sande et al. 2013 and Patel et al. 2013, these quiescent galaxies at z ^ 2 
grow inside-out, consistent with the expectations from minor mergers. 

4-3. Main sequence 

Recent studies have shown that there exists a correlation between star formation rate and 
stellar mass of galaxies at different redshifts {0< z <2), and it is usually called “main sequence” 
(Brinchmann et al. 2004; Daddi et al. 2007; Elbaz et al. 2007; Pannella et al. 2009; Rodighiero et 
al. 2011; Barro et al. 2013; Speagle et al. 2014; Whitaker et al. 2014). To investigate the essential 
star formation status of DRGs, we show the “main sequence” relation in the four redshift bins 
in Figure 8. In this hgure, the blue and red points represent star-forming and quiescent DRGs, 
classihed by the rest-frame UVJ colors, respectively. To specify the dependence of SFR on 
stellar mass, we also plot a linear-ht line (blue solid line) and its la dispersion (blue dashed 
line) for star-forming DRGs in each redshift panel. For comparison, we overplot the best-ht 
“main sequence” relation for z ^2 galaxies from Daddi et al. 2007 and z^l galaxies from Elbaz 
et al. 2007. From this hgure, we hnd a positive correlation between SFR and M* of DRGs in 
all four redshift bins. The slopes of the SFR—M* correlation of star-forming DRGs in diherent 
redshift panels are similar. At hxed stellar masses, the SFRs of galaxies on the SFR—M* 
correlation increase with redshifts, which indicates that star-forming DRGs were much more 
active on average in the past, due to greater abundance of gas in early universe. 

Figure 9 shows the dependence of specihc SFR (sSFR, sSFR = SFR/M^) on redshifts 
in four M* bins. The sSFR — Zph relation of star-forming and quiescent DRGs are shown in 
the left and right panel, respectively. In each panel, we divide DRGs into four M* bins (see 
Figure 9). For each M* bin, we calculate the (sSFR) (diamond symbols) and its standard 
deviation (error bars) in each redshift bin (Azph = 0.4). We hnd the sSFRs of DRGs increase 
with redshifts in all stellar mass bins, and DRGs with higher stellar masses generally have lower 
sSFRs. This is most likely due to the fact that the universe at early time is more active, and 
massive galaxies formed most of their stars earlier and slowed down more rapidly than their 
low-mass counterparts, star formation contributes more to the growth of low-mass galaxies 
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Fig. 8. SFR vs. M, diagram for DRGs in four redshift bins. Star-forming DRGs are plotted in 
blue points, quiescent DRGs are plotted in red points. The blue lines represent the best-fit of the main 
sequence of star-forming DRGs, and the blue dashed lines represent the Icr dispersion. For comparison, 
the best-fit lines from Daddi et al. 2007 and Elbaz et al. 2007 are also plotted. 
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Fig. 9. sSFR vs. Zph diagram for star-forming DRGs in left panel and quiescent DRGs in right panel. 
DRGs in different mass bins are plotted with different colors. The diamond symbols represent average 
sSFR value in redshift bin Azph = 0.4, and the error bars represent its Icr deviation. The histograms in 
the bottom panels show the redshift distribntions of DRGs. 

than high-mass galaxies. This result is in agreement with the “downsizing” mass-dependent 
quenching scenario (Kajisawa et ah 2010; Kajisawa et ah 2011). 

4.4- Star forming status and [3.6]-[24] color 

The total SFR of galaxy we used are determined via SFRtot=SFRuv,uncorr+SFRiR 
(Muzzin et ah 2013a). The SFRuv.uncorr comes from SED htting method and the SFRjr is 
derived by 24 /rm flux. Higher fraction of SFRjr in SFRtot means there are more UV photons 
absorbed by dust, and these energies are re-emitted from IR range. We plot the SFRiR/SFRtot 
vs. M* of all galaxies (left panel) and DRGs (right panel, red points for qDRGs, and blue points 
for sDRGs) in Figure 10. From this hgure we hnd there exists a trend in the SFRia/SFRtot 
vs. M* diagram of all galaxies: high-mass galaxies tend to have higher SFRi^/SFRtot ratios, 
indicating that massive galaxies tend to contain more dust and they are more metal rich, thus 
the extinction for UV lights are heavier than low-mass galaxies. Similar results are also found 
by Tadaki et al. 2013 in their Ha emitter (HAE) sample. We also hnd DRGs occupy the high- 
mass and high ratio of SFRia/SERtot region, which suggests the J — K color criterion effectively 
selects the massive dusty sources whose SFRir are dominant. As the SFRiR/SFRtot of most 
of DRGs great than 95%, it is suitable for DRGs using the 24 pm derived SFRjr to trace the 
whole star formation activity. 
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Fig. 10. SFRiR/SFRtot vs. stellar mass (M*) diagram for DRGs (right panel) and all galaxies in the 
COSMOS/UltraVISTA field (left panel). The dashed line corresponds to SFRiR/SFRtot = 95%. 


There are 29.4% quiescent and 49.8% star-forming all iG-selected galaxies in our sample 
have both [3.6] and [24] detection, and there are 37.3% quiescent and 87.7% star-forming DRGs 
have been detected at [3.6] and [24] bands simultaneously. The origins of the 24/im fluxes of 
the quiescent population can be due to processes unrelated to ongoing star formation, such as 
AGNs, cirrus dust heated by old stellar populations and circumstellar dust. Fumagalli et al. 
2014 treat each of these components separately and compare their contributions to Lir with the 
observed stacked values of Ljr, and conclude that various processes other than star formation 
can contribute to the observed mid-IR flux, thus the 24/im-derived SFRs for quiescent galaxies 
should be the upper limits. 

At the redshift z ~ 1 — 3, the IRAG 3.6 pm probes the rest-frame NIR range, which can 
be derived into stellar mass (Bell & de Jong 2001; Gole et al. 2001). The left panel of Figure 11 
shows the relationship of [3.6] vs. M* of all DRGs. We divide DRGs into four redshift bins, and 
hud there exists a strong correlation between [3.6] and M*. The linear-fit standard deviations 
of four redshift bins (from low to high) are 0.25, 0.20, 0.12 and 0.16, respectively. Similarly, the 
strong correlation between sSFRs and [3.6]-[24] colors in wide redshift range (1 < z < 3.5) are 
also found in the right panel of Figure 11. The standard deviations of linear £t of four redshift 
bins (from low to high) are 0.21, 0.19, 0.17 and 0.27, respectively. The average [3.6]-[24] color 
of quiescent all iF-selected galaxies is 0.78, while star-forming all iF-selected galaxies is 2.03. 
For our DRG sample, the average [3.6]-[24] color of quiescent DRG is 1.57, while star-forming 
DRG is 2.25. The [3.6]-[24] color of most of these 24pm-detected quiescent DRGs are generally 
consistent with but slightly redder than that of QGs presented by Brammer et al. 2009, Wang 
et al. 2012 and Huang et al. 2013, who identify QGs with the empirical [3.6]-[24] color criterion 
which is equivalent to a sSFR~ 10“^‘^(yr“^). Quiescent and star-forming populations separate 
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Fig. 11. Left panel: correlation between stellar mass and 3.6 magnitude of DRGs. Right panel: 
correlation between sSFR and [3.6]-[24] color of DRGs. DRGs in different redshift bins are plotted in 
different colors. Solid lines are best-fit lines for each redshift bin, dashed lines represent Icr deviation. 

obviously in both DRG and all iG-selected galaxy samples, which implies that [3.6]-[24] color 
is a good indicator for tracing star-forming status of galaxies. The similar results were also 
described in Huang et ah 2013. 

To demonstrate how the DRG selection picks up each type of galaxies, we show the frac¬ 
tion of quiescent galaxies and the fraction of dusty star-forming galaxies with SFRjR/SFRtot > 
0.95 selected by the DRG criterion as a function of redshift in Figure 12. Based on the stel¬ 
lar mass distribution of DRGs, we show the fractions for QGs and SFGs in three mass bins: 
10 < logio(M*/M 0 ) < 10.5, 10.5 < logio(M*/M 0 ) < 11 and 11 < logio(M*/M 0 ) < 12, which are 
represented by the three columns of Figure 12 from left to right. The panels in the top and 
bottom row are fractions for QGs and SFGs, respectively. From this hgure we hnd that (1) 
Higher fraction of dusty SFGs with SFRiR/SFRtot > 0.95 will be selected as DRG in higher 
mass bins, while the fractions of QGs remain nearly constant from low to high mass ranges; 
( 2 ) Both fractions of QGs and SFGs increase with redshift in each stellar mass bin, but the 
increase for QGs are more dramatic: very few QGs can be selected as DRG at z < 2 , while at 
z> 2 the fraction suddenly rises to a very high level. The result shown by this hgure indicates 
that heavy dust reddening oi J — K color exists among high-mass SFGs, while the J — K color 
of QGs is irrelevant to mass, because there is no dust extinction, the dramatic increase of the 
fraction of QGs is due to the sharp Balmer or 4000 A breaks of older stellar population shift 
into J-band at z> 2. 

4-5. Rest frame U-V color vs. stellar mass 

Galaxies present bi-modal distribution on the rest-frame color vs. mass diagram (GMD), 
and they show red rest-frame colors due to dust extinction or old stellar ages, thus the bi- 
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Fig. 12. Panel a), b) and c): Fraction of quiescent galaxies selected by the DRG criterion as a function of 
redshift. Panel d), e) and f): Fraction of dusty star-forming galaxies with SFRi^/SFRtot > 0.95 selected 
by the DRG criterion as a function of redshift. The fractions of QGs and SFGs are shown in three stellar 
mass bins (logio(M*/M0)=lO-lO.5, 10.5-11 and 11-12), and in each mass bin we divide them into six 
redshift bins (z=0.5-l, 1-1.5, 1.5-2, 2-2.5, 2.5-3 and 3-3.5) to calculate the fractions. 

modal distribution on CMD suggests two distinct types of galaxy formation histories (Bell et 
al. 2004; Cassata et al. 2008; Brammer et al. 2009). Talia et al. 2013 shows the good correlation 
between (f/ —5)rest and stellar mass, galaxies with redder {U — B)^est colors are more massive, 
no matter what kind of morphologies they are. Wang et al. 2012 hnd that disk galaxies tend 
to be more massive and have redder {U — ld)j,gg^ colors than irregular galaxies in their lERO 
sample. Diaz Tello et al. 2013 showed that with a particular (u — B) color, the high-redshift 
galaxies had on average higher sSFRs. 

Figure 13 shows the distribution of DRGs on {U — R)rest vs. M* diagram in six redshift 
bins. Star-forming and quiescent DRGs are plotted in blue and red, respectively. DRGs having 
visual morphological types introduced in Section 3.3 are also plotted in corresponding symbols. 
We hnd that the (f/ — R)rest color and M* of DRGs correlate very well in all redshift bins, which 
is in good agreement with Figure 10: more massive galaxies have redder {U — V)rest colors due 
to their larger internal dust extinctions or old stellar ages, which is also presented by Kajisawa 
& Yamada 2006. However, star-forming and quiescent DRGs cannot be separated efficiently 
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Fig. 13. Rest-frame {U — V) color vs. stellar mass diagram for DRGs in six redshift bins. Star-forming 
DRGs and quiescent DRGs are plotted in blue and red points, respectively. DRGs which have visual 
morphological types are also plotted in corresponding symbols. 

on CMD, they do not present clear bimodality, which is similar to what Wang et ah 2012 and 
Huang et ah 2013 have found. DRGs with similar stellar masses have bluer {U — V)iest colors at 
high redshifts, implying DRGs evolve with time, they gradually consume matter materials and 
eventually slow down their star forming activities (Kajisawa & Yamada 2005). We also hnd that 
the irregular DRGs dominate the reddest (f/ — R)rest color and high-mass end at high-redshifts, 
but gradually turn bluer and less massive than disk and spheroid DRGs toward low-redshifts. 
At z > 1.6, there are many blue spheroid DRGs, they are in the stage of compact star-forming 
galaxies, which are supposed to be the progenitors of compact QGs at high redshift. The 
hndings above may indicate a possible evolution scenario that massive disk-like DRGs merge 
with each other frequently by dissipative major merger at high redshifts, and produce massive 
dust-obscured irregular star-bursts, then through violent relaxation and angular momentum 
loss, they become spheroid-like DRGs and then quenched at later epoch. The hndings are also 
supported by Bedregal et ah 2013, Barro et ah 2013, Bassett et al. 2013 and Williams et ah 
2014. 
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Fig. 14. (a): Distribution of star-forming DRGs (blue) and quiescent DRGs (red) on ERO selection 
diagram. The selection criteria are introduced in Simpson et al. 2006 and Lane et al. 2007. (b): 

Distribution of DRGs on BzK selection diagram. The selection criteria are defined by Daddi et al. 2004. 
(c): Distribution of DRGs on lERO selection diagram. The selection criteria are defined by Yan et al. 
2004. (d): Distribution of DRGs on high-z ULIRG selection diagram. This selection criteria are defined 
by Huang et al. 2009. The insert panel show the distribution of DRGs on the Spitzer 24 pm magnitude. 
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5. Compare DRG to other color-selected high-redshift galaxies 

Over the past decade, many color selection techniques, which are designed based on 
rest-frame spectral energy distribution features, are applied to identify high-redshift galaxies 
(Steidel et ah 1996; Small et ah 1997; Franx et ah 2003; Daddi et ah 2004; Gronwall et al. 
2007; Yan et al. 2007; Dey et ah 2008; Dunne et al. 2009; Huang et al. 2009; Caputi et al. 
2012; Wang et al. 2012; Arcila-Osejo & Sawicki 2013; Ilbert et al. 2013). However, due to 
the small sizes of samples used in previous works, the overlaps between DRG and other color 
selected population have not been well studied. In this section, we use our large DRG sample 
to study the overlaps between DRGs and EROs, BzKs, lEROs and high redshift ULIRGs. 

First of all, we study the overlap between DRGs and EROs. To be consistent with Kong 
et al. 2006 and Lane et al. 2007, we use R — K color to select EROs from DRG sample. The color 
criterion (Rsubam — 77uKiDss)vega > 5.3 in Lane et al. 2007 can be converted to (i? —JF)ab > 3.65 
via i?AB = Rvega + 0.22 and TFab = ^Wega + 1.87. Figure 14 (a) shows the distribution of DRGs 
onto {R — K) vs. K diagram. We hnd that the distribution on K of quiescent DRGs is similar to 
star-forming DRGs. However, the distributions of quiescent and star-forming DRGs o\i R — K 
color are obviously different, quiescent DRGs have much redder R — K colors than star-forming 
DRGs, which attributes to the SED features of old stellar populations in quiescent DRGs. 69% 
DRGs can be selected as ERO, which represents the low redshift range of DRG population. 
Many DRGs especially sDRGs distribute at higher redshift are missed by the ERO selection. 

Secondly, we follow the procedure of McGracken et al. 2010 to perform B — z color 
calibration and correct the K band accounting for JFuitraViSTA — -^wmcam = 0.051 described 
in Muzzin et al. 2013a, to see how DRGs satisfy the BzK criteria. Figure 14 (b) shows the 
distribution of DRGs on the {z — K) vs. {B — z) plane. Star-forming DRGs are slightly bluer 
than quiescent DRGs on B — z color, but are obviously bluer on z — K color. 60.4% DRGs 
could be selected as sBzKs, 17.7% DRGs could be selected as pBzKs. Other DRGs mainly 
distribute in the redshift range oi z < 1.4, most (73%) of which are EROs. 

Thirdly, we analyze the overlap between DRGs and lEROs color criterion which is in¬ 
troduced by Yan et al. 2004. The purpose of lERO method is designed for selecting massive 
galaxies with mature stellar populations at z ^2 based on 4000 A Balmer break shifting beyond 
2 ;-band or steep slope of the rest-frame UV spectra of galaxies. Thus, a simple color cut 2;-[3.6] 
is designed to select high redshift galaxies. Young and blue galaxies cannot reach this threshold 
even at higher redshifts (Wang et al. 2012). Based on the simulation performed by continuous 
star formation (GSF) models in BG03 code, we obtain zacs — 2:subam = —0.11 to calibrate z- 
band. As shown in Figure 14 (c), distributions on [3.6]-[4.5] color of two population of DRGs 
are almost the same. However, star-forming DRGs are obviously bluer than quiescent DRGs on 
z-[3.6] color. There are 72% quiescent DRGs and 39% star-forming DRGs could be selected as 
lERO. Quiescent DRGs are more likely to be selected as lERO, the lERO color criterion rules 
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Fig. 15. Redshift distributions of DRGs that satisfy the ERO, BzK, lERO and high-z ULIRG color 
criteria. Red lines represent the distribution of DRGs satisfying the corresponding color criterion in each 
panel and black lines represent remanent DRGs. 

out a large proportion of star-forming DRGs, most of which are low-mass star-forming galaxies 
at low redshifts. 

Lastly, we perform comparison between DRG and high-z ULIRG color criteria introduced 
by Huang et al. 2009. This IRAG color selection method is designed for selecting massive, star¬ 
forming dominated ULIRGs at ^ ~ 2 based on the 1.6 pm stellar emission bump feature in 
galaxy SEDs. In 1.5 < 2 ; < 3, the four IRAG bands probe the rest-frame NIR bands where 
galaxy SEDs have similar shape, thus the IRAG colors are robust in determining redshift in 
this range (Huang et al. 2009). The color criterion is set as 0.05 < [3.6] — [4.5] < 0.4 and 
—0.7 < [3.6] — [8.0] < 0.5 based on the M82 SED model (Huang et al. 2004). Figure 14 (d) 
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Fig. 16. Distributions of other color-selected galaxies on the DRG selection diagram. Panel a), b), c), d) 
and e) shows the distribution of ERO, sBzK, pBzK, lERO and high-z ULIRG, respectively. 

shows the distribution of DRGs on the [3.6]-[8.0] vs. [3.6]-[4.5] diagram. We hnd that the 
peaks of DRGs on both two colors fall into the selection region designed by Huang et ah 2009, 
49.2% DRGs can be selected by this rectangle region, which reflects this selection method is 
efficient in determining redshifts. However, there are only 1.3% of DRGs satisfy their third 
cut, F 24 ^m > 0.5 mJy, which is correspond to 17.15 mag in AB, as shown in the insert panel 
of Figure 14 (d). Figure 15 shows the redshift distributions of DRGs, and the subsample of 
DRG satisfying other color criteria. We find that a large proportion of DRGs could also be 
selected as ERO, they mainly distribute at 2 ; ~ 1, most of which are dust reddened star-forming 
galaxies. The BzK method is very sharp in determining redshifts. Most of DRGs in our sample 
can be selected as sBzKs or pBzKs which distribute at z > 1.4, non-BzKs mainly distribute at 
z < 1.4. The lERO criterion introduced by Wang et ah 2012 and the high-z ULIRG criterion 
introduced by Huang et ah 2009 are all efficient in selecting 2 ; ~ 2 galaxies. 

Similarly, we show how other color-selected galaxies satisfy the DRG selection criterion 
in Figure 16. Panel a), b), c), d) and e) shows the distribution on the DRG selection J — K vs. 
K diagram for ERO, sBzK, pBzK, lERO and ULIRG, respectively. 24.3% EROs, 10.7% sBzKs, 
34.4% pBzKs, 64.1% lEROs and 15.8% high-z ULIRGs (not include the 24ynm flux criterion) 
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can be selected as DRG by the criterion J — K > 1.16. From the comparisons performed above, 
we can get a clear point of view in different color selection criteria: methods designed for 
selecting galaxies at z ~ 2 have many overlaps between each other. The DRG population is not 
a special galaxy group, the redshift distribution of them is broad, and most of them can also 
be selected by other color criterion. 

6. Summary 

In this paper, we have described the construction of a large DRG sample in the 
GOSMOS/UltraVISTA held. Based on multi-band photometry data, we analyze the physi¬ 
cal properties of DRGs. Our main conclusions are as follows. 

1. We select a sample of 4485 DRGs with {J — K)p^-q > 1.16 and JFab <23.4 mag, and 
identify 3725 star-forming DRGs and 760 quiescent DRGs based on rest-frame UVJ colors. We 
hnd the redshift distributions of DRGs are broad, ranging from 0.5 to 3.5 with a mean value 
of 1.79, but most of them (95.5%) distribute at 1 < z <?>. The stellar masses of DRGs mainly 
distribute at IO^^Mq — IO^^'^Mq. DRGs represent 70% of all galaxies in higher-mass bins at 
z > 2, indicating that the DRG selection criterion is sensitive to select massive galaxies at high 
redshift . 

2. We study the morphological properties of DRGs in our sample, employing data from 
HST WFG3/F125W and F160W imaging within the GOSMOS held. For DRGs distribute 
at 1 < z < 1.6 and 1.6 < 2 ; < 3, we measure Gini and M 20 from F125W and F160W images, 
respectively. The morphological results are consistent with our rest-frame UVJ color classih- 
cation: quiescent DRGs are generally compact and star-forming DRGs tend to have extended 
structures, implying that these galaxies have experienced diherent formation processes. 

3. The SFR and M* of star-forming DRGs present good “main sequence” relations in 
all redshift bins. Moreover, the sSFR of DRGs increase with redshift in all stellar mass bins 
and DRGs with higher stellar masses generally have lower sSFRs, which indicates that galaxies 
were much more active on average in the past, and star formation contributes more to the mass 
growth of low-mass galaxies than to high-mass galaxies, owning to massive galaxies formed 
most of their stars earlier and rapidly slow down than their low-mass counterparts. This result 
is in agreement with the “downsizing” scenario. 

4. There exists a trend in the SFRiR/SFRtot vs. M* diagram of all galaxies: high-mass 
galaxies tend to have higher SFRiR/SFRtot ratios, indicating that massive galaxies tend to 
contain more dust, thus the extinction for UV lights are heavier than low-mass galaxies. We 
hnd DRGs occupy the high-mass and high ratio of SFRjR/SFRtot region, which suggests the 
J — K color criterion effectively selects the massive dusty sources whose SFRir are dominant, 
and it is suitable for DRGs using [3.6]-[24] color to trace their sSFRs. The fraction of QGs 
and the fraction of SFGs with SFRiR/SFRtot > 0.95 selected by DRG criterion as a function of 
redshift indicate that the red J — K colors are mainly caused by dust reddening for SFGs and 
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the Balmer or 4000 A breaks shifting into the J-band at z > 2 for QGs. 

5. The ([/ — f4)rest color strongly correlates M* of DRGs in all redshift bins. Massive 
DRGs have redder (If — V)rest colors and the (17 — V)rest colors of DRGs become bluer at higher 
redshifts, suggesting high-mass galaxies have larger internal dust extinctions or older stellar 
ages and they evolve with time. We also hnd that the irregular DRGs dominate the reddest 
(U — V)rest color and high-mass end at high-redshifts, but gradually turn bluer and less massive 
than disk and spheroid DRGs toward low-redshifts. At z > 1.6, there are many blue spheroid 
DRGs. They are in the stage of compact star-forming galaxies which are supposed to be the 
progenitors of compact QGs at high redshift. These hndings may indicate a possible evolution 
scenario that massive disk-like DRGs merge with each other frequently, then through violent 
relaxation and angular momentum loss, they dually become spheroid-like galaxies. 

6. We study the overlaps between DRGs and EROs, BzKs, lEROs, high-z ULIRGs, hnd 
DRGs in our sample can be selected by other color criteria also. 69% DRGs can be selected as 
EROs, and most of them have low redshift. 60% DRGs can be selected as sBzKs, 18% DRGs 
can be selected as pBzKs, non-BzKs mainly distribute at z < 1.4, most (73%) of them are EROs. 
72% quiescent DRGs and 39% star-forming DRGs could be selected as lERO, quiescent DRGs 
are more likely to be selected as lERO, the lERO color criterion rules out a large proportion 
of star-forming DRGs, most of them are low-mass star-forming galaxies at low redshifts. 49% 
DRGs can be selected by the rectangle region of high- 2 ; ULIRG candidates and only 1.3% of 
DRGs in our sample satisfy the F 24 ^m > O.SmJy hux cut, which introduced by Huang et al. 2009. 
In turn, 24.3% EROs, 10.7% sBzKs, 34.4% pBzKs, 64.1% lEROs and 15.8% high-z ULIRGs 
(not include the 24pm hux criterion) can be selected as DRG by the criterion J — K > 1.16. 
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